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In this paper we present properties of relativistic and non-relativistic 
perfect hydrodynamical models. In particular we show illustrations of the 
fact that different initial conditions and equations of state can lead to the 
same hadronic final state. This means that alone from the hadronic ob- 
servables one cannot determine either of the above, one needs for example 
penetrating probes that inherit their properties from each timeslice of the 
evolution of the fireball. 

PACS numbers: 24.10.Nz 25.75.Ag 25.75.Ld 

1. Perfect fluid hydrodynamics 

Perfect fluid hydrodynamics is based on local conservation of entropy or 
density, energy and momentum, expressed by so-called conservation equa- 
tions. The fluid is perfect if the energy-momentum tensor is diagonal in 
the local rest frame, i.e. there are no shear stress, viscosity or heat conduc- 
tion effects. The conservation equations are closed by the equation of state, 
which gives the relationship between energy density e, pressure p. Typically 
e = up, where the proportionality "constant" k may depend on temperature 
T, which is connected to the density n and pressure p via p = riT in many 
solutions of the exactly solvable class, especially those discussed in following. 
In addition a bag constant B can be introduced with a modified equation 
of state of e — = K,{p + B), where i? = in the hadronic phase, non-zero 
in a deconfined phase. In this paper only the e = np case is discussed. 

Solutions describing non-relativistic flows that are applicable to relativis- 
tic heavy ion collisions are e.g. described in refs. [Il[2l[3], we will investigate 
a family of such solutions. 

Solving the relativistic equations is much harder. There are only a 
few exact solutions for these equations. One (and historically the flrst) 
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is the Landau-Khalatnikov implicit solution solution discovered more than 
50 years ago [H El [6]. This is a 1+1 dimensional solution, and has real- 
istic properties: it describes a 1+1 dimensional expansion, does not lack 
acceleration and predicts an approximately Gaussian rapidity distribution. 

Another renowned solution of relativistic hydrodynamics is the Hwa- 
Bjorken solution [7l[8l[9], which is a simple, 1+1 dimensional, explicit and 
exact, but accelerationless solution. This solution is boost-invariant in its 
original form, hence fails to describe the data [10(111]. However, the solution 
allowed Bjorken to obtain a simple estimate of the initial energy density 
reached in high energy reactions from final state hadronic observables. 

Important are solutions \12\ [T3] which are explicit and describe a rela- 
tivistic acceleration, i.e. combine the properties of the above solutions. We 
will investigate such a family of solutions. 



2. Investigated solutions 

We will investigate here an exact relativistic and an exact non-relativistic 
solution. We will use these solutions to extract information on the depen- 
dence of the final state on initial state, the parameters of the exact solution 
and the equation of state used. 



A non-relativistic hydro solution 

The below discussed non-relativistic solution describes a 3+1 dimen- 
sional ellipsoidally symmetric expansion [H [14], with an arbitrary k G R. 
The velocity field v is 

f X Y Z \ 
v= \j^r^,yry,-r, I , (1) 

where X, Y and Z are time-dependent principal axes of the expanding ellip- 
soid, X, Y and Z are their expansion rate versus time, and r = {rx,ry,rz) 
are the spatial coordinates. As for the thermodynamical quantities, one has 

XqYoZq ( rl rl rl \ 
n{r, t) = no^^ (^"^ " ^ " ^ j ' 

nr,t)=To(^Y\ (3) 



where Xq = X{to), Yq = Y{to), Zq = Z{to) are the principal axes at a 
given (arbitrarily chosen) time to, and uq = n(0, to)) ^0 = ^(0)*o)- This 
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represents a solution if the ordinary differential equations 

XX = YY = ZZ = a{^^^)"'^ (4) 

are fulfilled. Here a is an "acceleration parameter" , in the above mentioned 
solution a = TQ/m, but note, that with a generalization in the expression 
for n(r, t) a can be made arbitrary. 

A relativistic hydro solution 

Next being discussed is a family of analytic, explicit and simple so- 
lutions, which do not lack acceleration, and yield finite, realistic rapidity 
distributions [12]. In these solutions the velocity field is given by 

u'^ = 7(l,v) (5) 

with V being the d dimensional spherically symmetric velocity. The length 
of this vector is denoted by v, and this length can be expressed in the Rindler 
coordinates rj and r by 

V = tanh(Ar/) (6) 

where A is a kind of "acceleration parameter" , because relativistic accelera- 
tion u^d'^u'^ vanishes if A = 1, but is non- vanishing for a positive A different 
from 1. The pressure p and the temperature T are given by 

P = Po [—) ' (') 
T = To[^y'\ (8) 

The value of the constants A ( "acceleration parameter" ) , d (number of spa- 
tial dimensions) and k (adiabatic index of the fluid) are constrained, and 
different possible set of values yield different solutions. For example if A = 1, 
then d € R and k G II are arbitrary (this is corresponds to the Bjorken so- 
lution). In this case there is no acceleration. If however A = 2, i.e. there 
is substantial acceleration, d G R is arbitrary but k = d must be fulfilled. 
Another important case is, that if d = 1 and k = 1, then A G IR is arbitrary. 
This last case has a remarkably general velocity field: the A acceleration 
parameter can be arbitrary. On the other hand, this solution works only for 
d = 1 and k = 1, which is obviously a drawback. For a discussion of greater 
detail Refs. \12\ [T5] are recommended. 
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3. Dependence on the initial conditions and equation of state 

To calculate observables from the above solutions, freeze-out criteria 
have to be utilized. For example the simple condition of T(r, t) = Tfreeze-out 
can be used, or a condition where the freeze-out hypersurface is pseudo- 
orthogonal to the velocity field. This is an important aspect of hydrody- 
namics, that a freeze-out condition has to be chosen, but here it does not 
need to be discussed: we fix the hydrodynamical final state and investigate 
only those hydrodynamical evolutions that lead exactly to the same final 
state (i.e. to the same hadronic observables). 

First, let us see relativistic solutions of section [2l In Figs. [Hand [2] the 
spatial temperature distribution is shown at difi'erent times, for different ini- 
tial conditions and difi'erent A and/or k, values. All the plotted solutions go 
to the same final temperature at mid-rapidity [rj = 0), although their initial 
condition, acceleration parameter and/or equation of state is diff^erent. 

Then let us see non-relativistic solutions of section [2l In Fig. [3] time 
evolution of the principal axes of the expanding ellipsoid, their expansion 
rates and the temperature is shown, for difi'erent equations of state and 
different initial conditions. Here parameters were chosen so that the final 
states of all cases are the same: Tq = 200 MeV, X = 11 fm, y = 12 fm, 
Z = 13 fm, Xo = 0.6, io = 0.5, Zq = 0.7. 

Conclusions 

We have seen that the same hydrodynamical final state can be achieved 
with different solutions, equations of state or initial conditions A very im- 
portant consequence is also, that the hydrodynamic scaling laws observed 
in relativistic heavy ion collisions [16] do not depend on the initial con- 
ditions or the equation of state separately - one needs to restrict one in 
order to extract information on the other. These examples show that in or- 
der to extract the initial conditions or the equation of state from hadronic 
(final-state) observables, one needs further restrictions or experimental con- 
straints. With penetrating probes one has access to the earlier times of the 
evolution thus can restrict results on initial conditions or equation of state. 
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Fig. 1. These plots show the temperature distribution of a c? = 1 exact relativistic 
solutions from section[21 for different A parameters and different initial temperature 
distributions. All approximate the same final temperature distribution {T « Tq). 
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Fig. 2. Temperature distribution is shown for a d = 1 exact relativistic solution 
from section [5] with different initial temperatures and equations of state (k). All 
approximate the same final temperature distribution. 
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Fig. 3. Time evolution of the principal axes of the expanding ellipsoid {X,Y,Z), 
their expansion rates {X, Y, Z, denoted by Xv etc.), as well as that of the temper- 
ature (T) is plotted, for different equations of state (k). It is evident that different 
equations of states combined with different initial conditions can lead to the same 
final state. The solution used here is the exact non-relativistic solution of section[2l 



